arXiv:1508.06562v2 [hep-ex] 9 Jun 2016 


Energy scan of the e+e —)■ hb{nP)n~^7T (n = 1,2) cross sections and evidence for 

T(11020) decays into charged bottomoninm-like states 

1. Adachi,^^’® H. Aihara,®® D. M. Asner,^^ H. Atmacan,^® V. Aulchenko,®’®° T. Aushev,"^^ R. Ayad,®^ 

1. Badhrees,®^’A. M. Bakich,®° E. Barberio,®^ P. Behera,^^ V. Bhardwaj,®® B. Bhuyan,^® J. Biswal,^^ 

A. Bobrov,®’®® A. Bondar,®’®® G. Bonvicini,^® A. Bozek,^^ M. Bracko,®®’^^ T. E. Browder,®^ D. Cervenkov,^ 

V. Chekelian,®® A. Chen,^^ B. G. Cheon,®® K. Chilikin,®®’^® R. Ghistov,®®’-®® V. Ghobanova,®® S.-K. Ghoi,®® 

Y. Choi,®® D. Cinabro,®"® J. Dalseno,®®’®® M. Danilov,^®’®® N. Dash,®® Z. Dolezal,^ A. Drutskoy,®®’^® 

S. Eidelman,®’®® D. Epifanov,®® T. Ferber,® B. G. Fulsom,®® V. Gaur,®® A. Garmash,®’®® R. Gillard,®"® Y. M. Goh,®® 
P. Goldenzweig,®^ B. Golob,®®’®® D. Greenwald,®"® T. Kara,®®’® K. Hayasaka,®”® H. Hayashii,^® W.-S. Hou,^® 
C.-L. Hsu,®®’ K. Inami,'®® G. Inguglia,® A. Ishikawa,®® Y. Iwasaki,®® I. Jaegle,®® T. Julius,®®’ K. H. Kang,®® 

P. Katrenko,"®®’®® D. Y. Kim,®®’ H. J. Kim,®® J. B. Kim,®®" K. T. Kim,®®’ M. J. Kim,®® S. H. Kim,®® Y. J. Kim,®® 

K. Kinoshita,® P. Kodys,'® S. Korpar,®®’®® D. Kotchetkov,®® P. Krokovny,®’®® T. Kuhr,®® A. Kuzmin,®’®® 

Y.-J. Kwon,®’® J. S. Lange,® C. H. Li,®®’ H. Li,®® L. Li,®^ L. Li Gioi,®® J. Libby,®®’ D. Liventsev,'^®’®® M. Lubej,®® 

T. Luo,®® M. Masuda,®®’ T. Matsuda,®® D. Matvienko,®’®® K. Miyabayashi,^® H. Miyata,^® R. Mizuk,®®’"®®’-®! 

G. B. Mohanty,®® A. Moll,®®’®® E. Nakano,®® M. Nakao,®®’® T. Nanut,®® K. J. Nath,®® K. Negishi,®® M. Niiyama,®® 

N. K. Nisar,®®’® S. Nishida,®®’® S. Ogawa,®® S. Okuno,®® S. L. Olsen,®® Y. Onuki,®® P. Pakhlov,®®’^® 

G. Pakhlova,®®’^® B. Pal,® C. W. Park,®® H. Park,®® S. Paul,®^ T. K. Pedlar,®'® R. Pestotnik,®® M. Petrie,®® 

L. E. Piilonen,"^® C. Pulvermacher,®® M. Ritter,®® Y. Sakai,®®’® S. Sandilya,® T. Sanuki,®® V. Savinov,®® 

T. Schliiter,®® O. Schneider,®® G. Schnell,®’®'® C. Schwanda,®® Y. Seine,'®® D. Semmler,® K. Senyo,"^® 

O. Seon,®® M. E. Sevior,®®’ V. Shebalin,®’®o T.-A. Shibata,®® J.-G. Shiu,^® B. Shwartz,®’®® F. Simon,®®’®® 

E. Solovieva,®®’®® M. Staric,®® J. Stypula,®"^ T. Sumiyoshi,"^® M. Takizawa,®® U. Tamponi,®®’’^® K. Tanida,®® 

Y. Teramoto,®® 1. Tikhomirov,K. Trabelsi,®®’® M. Uchida,®® T. Uglov,®®’4® Y. Unno,®® S. Uno,®®’® 

P. Urquijo,®®’ Y. Usov,®’®® C. Van Hulse,® G. Varner,®® V. Vorobyev,®’®® C. H. Wang,^® M.-Z. Wang,'®® 

P. Wang,®® X. L. Wang,"^® Y. Watanabe,®® K. M. Williams,’^® E. Won,®"^ J. Yamaoka,®® Y. Yamashita,®® 

J. Yelton,"^ C. Z. Yuan,®® Z. P. Zhang,®^ V. Zhilich,®’®® V. Zhukova,®® V. Zhulanov,®’®® and A. Zupanc®®’®® 

(The Belle Collaboration) 

^Aligarh Muslim University, Aligarh 202002 
® University of the Basgue Country UPV/EHU, 48080 Bilbao 
^Budker Institute of Nuclear Physics SB RAS, Novosibirsk 630090 
* Faculty of Mathematics and Physics, Charles University, 121 16 Prague 
^University of Cincinnati, Cincinnati, Ohio 45221 
^Deutsches Elektronen-Synchrotron, 22607 Hamburg 
"^University of Florida, Cainesville, Florida 32611 
^Justus-Liebig-Universitdt Ciefien, 35392 Giefien 
^SOKENDAI (The Graduate University for Advanced Studies), Hayama 240-0193 
Gyeongsang National University, Chinju 660-701 
Hanyang University, Seoul 133-791 
®® University of Hawaii, Honolulu, Hawaii 96822 
^^High Energy Accelerator Research Organization (KEK), Tsukuba 305-0801 
IKERBASQUE, Basgue Foundation for Science, 48013 Bilbao 
Indian Institute of Technology Bhubaneswar, Satya Nagar 751007 
Indian Institute of Technology Guwahati, Assam 781039 
Indian Institute of Technology Madras, Chennai 600036 
Indiana University, Bloomington, Indiana 47408 
Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049 
Institute of High Energy Physics, Vienna 1050 
^^INFN - Sezione di Torino, 10125 Torino 
Stefan Institute, 1000 Ljubljana 
Kanagawa University, Yokohama 221-8686 

‘^^‘Institut fiir Experimentelle Kernphysik, Karlsruher Institut fur Technologic, 76131 Karlsruhe 
King Abdulaziz City for Science and Technology, Riyadh 1144^ 

Korea Institute of Science and Technology Information, Daejeon 305-806 
Korea University, Seoul 136-713 
Kyoto University, Kyoto 606-8502 
Kyungpook National University, Daegu 702-701 
Ecole Polytechnigue Federate de Lausanne (EPFL), Lausanne 1015 


2 


^^P.N. Lebedev Physical Institute of the Russian Academy of Sciences, Moscow 119991 
Faculty of Mathematics and Physics, University of Ljubljana, 1000 Ljubljana 
Ludwig Maximilians University, 80539 Munich 
Luther College, Decorah, Iowa 52101 
University of Maribor, 2000 Maribor 
Max-Planck-Institut fur Physik, 80805 Miinchen 
School of Physics, University of Melbourne, Victoria 3010 
Middle East Technical University, 06531 Ankara 
University of Miyazaki, Miyazaki 889-2192 
Moscow Physical Engineering Institute, Moscow 115409 
Moscow Institute of Physics and Technology, Moscow Region IfllOO 
Graduate School of Science, Nagoya University, Nagoya 464-8602 
'^^Nara Women’s University, Nara 630-8506 
National Central University, Chung-li 32054 
National United University, Miao Li 36003 
Department of Physics, National Taiwan University, Taipei 10617 
H. Niewodniczanski Institute of Nuclear Physics, Krakow 31-342 
^^Nippon Dental University, Niigata 951-8580 
Niigata University, Niigata 950-2181 
^^Novosibirsk State University, Novosibirsk 630090 
Osaka City University, Osaka 558-8585 
Pacific Northwest National Laboratory, Richland, Washington 99352 
University of Pittsburgh, Pittsburgh, Pennsylvania 15260 
University of Science and Technology of China, Hefei 230026 
Seoul National University, Seoul 151-742 
Show a Pharmaceutical University, Tokyo 194-8543 
Soongsil University, Seoul 156-743 
University of South Carolina, Columbia, South Carolina 29208 
Sungkyunkwan University, Suwon 440-746 
School of Physics, University of Sydney, New South Wales 2006 
Department of Physics, Paculty of Science, University of Tabuk, Tabuk 71451 
^^Tata Institute of Pundamental Research, Mumbai 400005 
Excellence Cluster Universe, Technische Universitdt Miinchen, 85748 Garching 
Department of Physics, Technische Universitdt Miinchen, 85748 Garching 
Toho University, Funabashi 274-8510 
Department of Physics, Tohoku University, Sendai 980-8578 
Earthquake Research Institute, University of Tokyo, Tokyo 113-0032 
Department of Physics, University of Tokyo, Tokyo 113-0033 
Tokyo Institute of Technology, Tokyo 152-8550 
Tokyo Metropolitan University, Tokyo 192-0397 
University of Torino, 10124 Torino 

Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 
Wayne State University, Detroit, Michigan 48202 
Yamagata University, Yamagata 990-8560 
Yonsei University, Seoul 120-749 

Kobayashi-Maskawa Institute, Naqoya University, Nagoya 464-8602 
(Dated: August 11, 2016) 

Using data collected with the Belle detector at the KEKB asymmetric-energy e''"e“ collider, we mea¬ 
sure the energy dependence of the e''"e“ —>■ hb{nP)n^'K~ {n = 1, 2) cross sections from thresholds 
up to 11.02 GeV. We find clear T(10860) and T(11020) peaks with little or no continuum contri¬ 
bution. We study the resonant substructure of the T(11020) —>■ hb{nP)'K'^'K~ transitions and find 
evidence that they proceed entirely via the intermediate isovector states .^6(10610) and ^^(lOBSO). 
The relative fraction of these states is loosely constrained by the current data: the hypothesis that 
only Zi,(10610) is produced is excluded at the level of 3.3 standard deviations, while the hypothesis 
that only Zb{10650) is produced is not excluded at a significant level. 

PACS numbers: 14.40.Rt, 14.40.Pq, 13.66.Bc 


Heavy quarkonium is a bound state of the cc or bb 
quarks. The heavy quarks are moving relatively slowly; 
therefore, a non-relativistic approximation based on the 
interaction potential accurately describes the basic prop¬ 


erties of this system [l[. The first state that did not fit 
potential model expectations was observed in 2003 by 
Belle Q ; since then, almost twenty such states have been 
reported Q. They correspond to high excitations and 
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have masses above the DD or BB thresholds. 


Many quarkonium-like states were found in the energy 
scans of the cross sections of e+e” annhilation into con¬ 
ventional quarkonia and light hadrons. Among these are 
the y(4008) and F(4260) in J/'(/'7r+7r“ Q, the F(4360) 
and Y (4660) in ii){2S) tt+tt” Q, the ■!/)(4040) and ^(4160) 
in and possibly the F(4220) in 0]. The 

partial widths of the corresponding transitions are much 
higher than expected for conventional quarkonia Q . Sur¬ 
prisingly, the peaks observed in the cross sections de¬ 
pend on the final states. In other words, each such 
charmonium-like state decays to only one channel with 
charmonium. To explain this “selectivity”, a hadrochar- 
monium notion is introduced : a bound state of a char¬ 
monium and a light hadron. Such a system decays pre¬ 
dominantly into its constituents. 

Recent energy scans of the e+e“ —T(nS')7r+7r“ 
(n = 1,2,3) cross sections by Belle EEl show that 
the situation is different in the sector of bottomonium- 
like states: all of the cross sections exhibit peaks of the 
T(10860) and T(11020) resonances that are also seen in 
the total hadronic cross section. The observed decay 
patterns of T(10860) and T(11020) agree with the ex- 
pectations for a mixture of the B^^^B^^^ molecule and 
conventional bottomonium: open flavor channels domi¬ 
nate, while channels with quarkonium have anomalously 
high partial widths [l^. The striking difference between 
charmonium-like and bottomonium-like states is not yet 
understood. Further scans in the bottomonium region 
are therefore of high importance. In this Letter, we re¬ 
port the first energy scan of the e+e“ —)■ /ib(nP)7r+7r“ 
(n = 1,2) cross sections. We find clear T(10860) and 
T(11020) peaks without a significant continuum contri¬ 
bution. 


To date, the e+e“ —>■ hb{nP)Tr~^Tr~ processes were seen 
only at a single energy near the T(10860) peak [l^ . 
They were found to proceed entirely via the intermediate 
isovector states Zb(10610) and Zt,(10650) that are situ¬ 
ated near the BB* and B*B* thresholds fl3 a nd likely 
have corresponding molecular structures [l4|. Here, 
we report on the resonant substructure study of the 
T(11020) —>• /ib(nP)7r+7r“ decays, where we find first 
evidence for intermediate states. Hereinafter, the 
T(10860) and T(11020) are referred to, for brevity, as 
the T(55') and T(6S') according to the potential model 
assignment. 


We use 121.4 fb“^ of on-resonance T{5S) data taken 
at three energies close to 10.866 GeV, as well as 1 fb“^ 
of data taken at each of 19 different energies between 
10.77 and 11.02 GeV. These data were collected with 
the Belle detector 0 at the KEKB asymmetric-energy 
e+e“ collider [17 1. 


The e“'"e —^ /ib(nP)7r+7r processes are recon¬ 
structed inclusively based on the tt+tt” missing mass, 
A7miss(7r7r) = \/(Ac.m. - E*^Y -Pll, where Ac.m. is the 


center-of-mass (c.m.) energy and E*^ and are the 
energy and momentum of the 7r“'"7r“ pair as measured in 
the c.m. frame. The c.m. energy is calibrated using the 
e+e" —?► T(n5)7r+7r“ —>■ /r+^“7r+7r“ and e“''e“ —?► 
processes, as described in Ref. [llj. This analysis closely 
follows previous Belle publications 1^ 14, Iq. 

We use a general hadronic event selection with re¬ 
quirements on the position of the primary vertex, track 
multiplicit y, and the total energy and momentum of 
the event [l9|. These criteria suppress Bhabha, pT^ 
r+r“, two-photon and beam-gas processes. Gontinuum 
e+e” qq {q = U, d, s, c) events have jet-like shapes in 
contrast to the spherically symmetric signal events and 
are suppressed by a requirement on the ratio of the sec¬ 
ond to zeroth Fox-Wolfram moments: i ?2 < 0.3 [l^. We 
only consider positively identified 7r“''7r“ candidates that 
originate from the interaction point region. 

The measurements of the cross sections are performed 
with an additional requirement on the single-pion 
missing mass, 

10.59 GeV/c^ < < 10.67GeV/c^ (1) 


which selects signal events proceeding via the interme¬ 
diate Zb{10610) or Zb(10650) states. We combine the 
A7miss(7r7r) distribution for tt’*' satisfying © and that for 
7T~ satisfying O- The 7r’''7r pairs with both Mmiss(7r~'") 
and Mniiss(7r“) in the Zj, mass window are counted twice. 
(If they are counted only once, the combinatorial back¬ 
ground develops a dip slightly above the hh{2P) signal, 
making the background parameterization difficult.) We 
take the double entries into account by correcting the 
errors of the Mmiss(7r7r) histogram and, based on Monte 
Garlo (MG) simulation, the hb{2P) signal yields. 

We fit the Mmiss(7r7r) distribution in the hb{lP) and 
hb{2P) intervals, defined as 9.8GeV/c^ — lO.OGeV/c^ 
and 10.17GeV/c^ — 10.34 GeV/c^, respectively. The 
fit function is the sum of the hb(nP) signal and 
combinatorial- and peaking-background components. 
The shapes of the hb{nP) signals are determined by con¬ 
volving the probability density of the initial state radia¬ 
tion (ISR) process with the experimental resolution, de¬ 
scribed by a Gaussian. We use the ISR probability, cal¬ 
culated up to the second order 2l| , and take into account 
the energy dependence of the e+e“ ^ hb{nP)n'^Tr~ cross 
sections using an iterative procedure. The resolution is 
determined using the exclusively reconstructed decays 
T(5iS') T(nS')7r+7r“, T{nS) to be (6.84 ± 

0.13) MeV/c^ for the hb{lP) and (6.15±0.22) MeV/c^ for 
the hb{2P). The resolution is dominated by c.m. energy 
smearing. The hb(nP) masses are fixed at the previous 
Belle measurement [l8|. We normalize the signal den¬ 
sity functions in such a way that the measured hbinP) 
yields include the ISR correction, 1 + i5isr, and can be 
used directly to measure the Born cross sections. The 
combinatorial background is described by a fourth-order 
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Chebyshev polynomial in both fit intervals. The order is 
chosen by maximizing the confidence level of the fit. 

Using MC simulation, we find that combining a ran¬ 
dom pion that satisfies the Z}, mass requirement and a 
signal pion from —>■ hi,(nP)TT produces a broad bump 
under the hb{nP) signal. This background is incorpo¬ 
rated within the combinatorial background and results 
in minor corrections in the hb{lP) and hb{2P) yields of 
0.99 ± 0.01 and 0.995 ± 0.005, respectively. The tt+tt” 
pairs originating from the T{2S) —>■ T(lS')7r+7r“ transi¬ 
tions with the T{2S) produced inclusively or via ISR re¬ 
sult in a peak at Ec.m. — [wx( 2 S) ~ nT-r^is)] that is inside 
the hb{2P) fit interval for the c.m. energies close to the 
T(55'). The shape of this peaking background is found 
to be a Gaussian with cr = 11 MeV/c^. Its normalization 
is floated in the fit. 

To determine the reconstruction efficiency, we use 
phase-space-generated MC, weighted in Mi„iss(7r) accord¬ 
ing to the fit results for the T(55') —>■ hb(lP)7r+7r“ tran¬ 
sitions 0 and in angular variables according to the ex¬ 
pectations for the Zb spin-parity = 1+ [i^. The 
efficiencies for the hb(lP)7r+7r“ and /ib(2P)7r+7r“ chan¬ 
nels are in the range 40 — 55% and 35 — 50%, respectively; 
they rise with c.m. energy. At the lowest energy point, 
there is a drop of efficiency by a factor of two since this 
point is close to the kinematic boundary and the pion 
momenta are low. 

At each energy, the Born cross section is determined 
according to the formula: 

cr^(e+e“ /if,(nP)7r+7r-) = (2) 

where N is the number of signal events determined from 
the Mmiss(7r7r) fit that includes the ISR correction, L 
is the integrated luminosity, e is the reconstruction ef¬ 
ficiency and |1 — np is the vacuum polarization correc¬ 
tion [^, which is in the range 0.927 — 0.930. The re¬ 
sulting cross sections are shown in Fig. [TJ The cross sec¬ 
tions, averaged over the three high statistics on-resonance 
points at Pc.m. = (10865.6 ± 2.0) MeV, are 

cr®(e+e“ ^ /ib(lP)7r+7r-) = 1.66 ± 0.09 ± 0.10 pb, (3) 

cr^(e+e- hb{2P)TT+TT-) = 2.70 ± 0.17 ± 0.19pb. (4) 

The ratio of the cross sections is 0.616 ± 0.052 ± 0.017. 
Here and elsewhere in this Letter, the first uncertainties 
are statistical and the second are systematic. 

The systematic uncertainties in the signal yields origi¬ 
nate from the signal and background shapes. The relative 
uncertainty due to the Mniiss(7’‘7’‘) resolution is correlated 
among different energy points and is equal to 1.4% for 
the hb{lP) and 3.3% for the hb{2P). The uncertainties 
due to the hb{nP) masses and ISR tail shapes are found 
to be negligible. To estimate the background-shape con¬ 
tribution, we vary the fit interval limits by about 50 MeV 



FIG. 1. (colored online) The cross sections for the e''"e“ 
hb{lP)'K^TT~ (top) and e^e~ —>■ hb{2P)TT^Tv~ (bottom) as 
functions of c.m. energy. Points with error bars are the data; 
outer error bars indicate statistical uncertainties and inner 
red error bars indicate uncorrelated systematic uncertainties. 
The solid curves are the fit results. 


and the polynomial order for each fit interval. The cor¬ 
responding uncertainties are considered uncorrelated and 
are 1.1% and 2.5% for the on-resonance cross sections in 
Eqs. ([3]) and ([4]), respectively. 

A relative uncertainty in the efficiency contributes to 
the correlated systematic uncertainty. An uncertainty 
due to the Zb mass requirement of i!ii;g% is estimated 
by varying the Zb parameters by Tier and taking into 
account correlations among different parameters. The 
efficiency of the R 2 requirement is studied using inclu¬ 
sively reconstructed T(5S') ^ T(nS')7r+7r“ decays. We 
find good agreement between data and MC and assign 
the 5% statistical uncertainty in data as a systematic un¬ 
certainty due to the R 2 requirement. Finally, we assign 
a 1% uncertainty per track due to possible differences in 
the reconstruction efficiency between data and MC. 

An uncertainty in the luminosity of 1.4% is primar¬ 
ily due to the simulation of Bhabha scattering that is 
used for its determination and is correlated among energy 
points. We add in quadrature all the contributions to find 
the total systematic uncertainties shown in Eqs. m and 
O- The values of the cross sections for all energy points 
are provided in Ref. [ 2 ^. 

The shapes of the hh(lP)7r+7r“ and hb{2P)Tr~^Tr~ cross 
sections look very similar. They show clear T(55') and 
T(6S') peaks without significant continuum contribu¬ 
tions. We perform a simultaneous fit of the shapes, 
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adding in quadrature the statistical and uncorrelated sys¬ 
tematic uncertainties at each energy point. We use the 
coherent sum of two Breit-Wigner amplitudes: 

A. $n(s) |FBw(s,M5,r5) + ae*-^FBw(s,M6,r6)P, (5) 

where s = , ‘hra(s) is the phase space calculated 

numerically, taking into account the measured Zt line 
shape [3, and Fbw{s,M,T) = MT/{s - IVP + iMT) 
is a Breit-Wigner amplitude. The fit parameters M 5 , 
Ts, Me, Tg, a and (j) are common for the two channels, 
while only the normalization coefficients An are differ¬ 
ent. Equation (O is convolved with the E^.m. resolution 
of (5.0 ± 0.4) MeV, which is found using exclusively re¬ 
constructed T(55') —> T(nS') 7 r+ 7 r“ events. The fitted 
functions are shown in Fig. [T] The confidence level of 


the fit is 93%. The fit results are: 

Ms = (10884.7t^ MeV/c^ (6) 

= (7) 

Mg = (10999.0t^j;^^;^) MeV/c^ (8) 

rg = (27;27+^5)MeV, (9) 


a = and </> = (0.lti(;8 ± 0.3) tt. (10) 


The measured masses and widths agree with the results 
of the T(n 5 ') 7 r+ 7 r“ scan 111 . 

The first error in the fit results is not purely statisti¬ 
cal but includes uncorrelated systematic uncertainties in 
the cross sections. The contributions of other considered 
sources are listed in Table ID 


TABLE I. The systematic uncertainties in the T(5S') and 
T(6S) masses (in MeV/c^), widths (in MeV), amplitude a, 
and phase 0 (in units of tt). 



Ms 

Fs 

Me 

Fe 

a 

0 

Fit model 

+8.9 

-0.1 

+ 0.4 
-19.1 

+ 16.7 
- 0.0 

+ 0.0 
-11.5 

+0.12 

-0.00 

+0.09 

-0.00 

Zb substructure 

+0.2 

-0.0 

+0.0 

-0.2 

+0.1 

-0.0 

+0.7 

-0.0 

+0.11 

-0.00 

+0.00 

-0.29 

^/s scale 

1.0 

1.0 

+3.0 

-1.0 

+4.7 

-1.0 

+0.00 

-0.10 

+0.25 

-0.00 

Resolution 

0.0 

+0.3 

-0.2 

0.1 

0.6 

0.0 

+0.01 

-0.00 

Total 

+8.9 

+ 1.1 

+ 16.9 

+ 4.8 

+0.17 

+0.27 

-1.0 

-19.1 

- 1.0 

-11.5 

-0.10 

-0.29 


To study systematic uncertainties due to the fit model, 
we introduce a non-resonant continuum amplitude, 

The significance of this contribution is only l.bcr. How¬ 
ever, the shifts in the fit results are large, and this is 
the dominant source of systematic uncertainty. We also 
consider the possibility that the parameters a and cj) are 
different in the /it,(lP) 7 r+ 7 r“ and hg( 2 P) 7 r+ 7 r“ channels. 
We find that the values in the two channels agree and 
the shifts in masses and widths are small. Using MC 
pseudo-experiments, we find that there is no significant 
fit bias. 

If the resonant substructures of the T(55') and T(65') 
decays are different, the T(55') and T(65') amplitudes 


in ((S)) are not fully coherent, and the interference term 
is suppressed by a decoherence factor k [ll|. If only 
Zb(10610) is produced at the T(65'), k is calculated nu¬ 
merically to be 0.62; if only Zb(10650) is produced, k is 
0.80. We introduce these factors in the fit and take into 
account that the efficiency of the Zb mass requirement is 
smaller for a single Zb state compared to two Zb states 
by 12 % since the two Zb states interfere destructively 
outside their signal region. 

We account for an uncertainty in the Pc.m. scale and 
the uncertainty in the Pc.m. resolution. We add in 
quadrature the contributions of the various sources to 
determine the total systematic uncertainties. 

To study the resonant substructure of the T( 6 S') —>■ 
hb{nP)'K^'K~ transitions, we combine the data samples 
of the five highest-energy points. The fits to the corre¬ 
sponding Mmiss( 7 r 7 r) spectra in the hb{lP) and hb{2P) 
intervals are shown in Figs. [2] and El The hb{nP) sig- 



9.8 9.85 9.9 9.95 10 

GeV/c^ 


FIG. 2. (colored online) The Mmiss(7r7r) spectrum in the 
hh{lP) region for the combined data samples of five energy 
points near the T(6S). In (a) the data are the points with er¬ 
ror bars with the fit function (solid curve) and background 
(red dashed curve) overlaid. (b) shows the background- 
subtracted data (points with error bars) with the signal com¬ 
ponent of the fit overlaid (solid curve). 


nal density functions are determined by averaging over 
the data samples that are combined; we use weights pro¬ 
portional to the integrated luminosity and the cross sec¬ 
tion at each energy. The confidence levels of the fits are 
50% and 52%, respectively. From Wilks’ theorem [3 - 
we find that the significances of the hb{lP) and hb{2P) 
signals are 3.5cr and 5.3(T, respectively, including sys¬ 
tematic uncertainty, determined by varying the poly¬ 
nomial order. Thus, we find the first evidence for the 
T( 6 S') —I- hb{lP)Tr~^TT~ transition and observe for the first 
time the T(65) hb( 2 P) 7 r+ 7 r“ transition. 

We release the requirement of an intermediate Zb and 
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10.2 ^ ^ 1o!25 ^ ^ 10.3 


M^,3(7iV), GeV/c^ 

FIG. 3. (colored online) The MmissiTm) spectrum in the 
hb{2P) interval for the combined data samples of five energy 
points near the T(65). The legend is the same as in Fig. (2] 


fit the Mi„iss(7''7r) spectra in bins of Mmiss(7r) to measure 
the hi,{nP)TT^TT~ yields as functions of Mmiss(7i‘)- The 
distribution of the phase-space-generated signal events 
in the Mmiss(7r“'") vs. Miniss(7r“) plane has the shape 
of a narrow slanted band; each structure at high val¬ 
ues of Mi„iss(7r^) produces a “reflection” at small values 
of We combine the Miniss(7r7r) spectra for 

the corresponding Mmiss (tJ'"'’) and bins and 

consider the upper half of the available Mi„iss(7r) range. 
Thereby, we consider all signal events and avoid dou¬ 
ble counting. The yields, corrected for the reconstruc¬ 
tion efficiencies, are shown in Fig. |4l The data are not 




FIG. 4. The efficiency-corrected yields of hb{lP)'K'^TT~ (a) 
and hb{2P)n'^'K~ (b) as functions of Miniss(7r) for the com¬ 
bined data samples of five energy points in the T {6S) region. 
Points represent data; the solid histogram represents the fit 
result with the Zb signal shape fixed from the T (SS") analysis; 
the dashed histogram represents the result of the fit with a 
phase space distribution. 


distributed uniformly in phase space; they populate the 


Zh(10610) and Zh(10650) mass region. We fit the data 
to a shape where the ^{,(10610) and ^{,(10650) parame¬ 
ters are fixed to the T(5S') —>■ Zbir —>■ /ib(lP)7r+7r“ re¬ 
sult and the non-resonant contribution is set to zero 141. 
Such a model describes the data well: the confidence 
levels of the fits are 65% and 77% for the hb{lP) and 
hb{2P), respectively. The phase space hypothesis is ex¬ 
cluded relative to this model at the 3.6a and A.5a levels 
in the /ib(lP)7r+7r“ and hb(2P)7r+7r“ channels, respec¬ 
tively. The single Z'h(10610) hypothesis is excluded at 
the 3.3cr level in the ft,b(lP)7r+7r“ channel, while the sin¬ 
gle Zb(10650) hypothesis cannot be excluded at a signif¬ 
icant level. In the /i6(2P)7r+7r“ channel, the Zb(10610)^ 
and Zb(10650)^ signals overlap with the Zt,(10650)^ and 
Zb(10610)^ reflections, respectively, which obscures the 
determination of the relative yields. The exclusion levels 
are determined using pseudo-experiments from the 
differences of the two hypotheses being compared, and 
include systematic uncertainty. 

In conclusion, we have measured the energy depen¬ 
dence of the e+e“ —/if,(nP)7r+7r“ (n = 1,2) cross sec¬ 
tions. We find two peaks corresponding to the T(5S') and 
T(65') states and measure their parameters, which agree 
with the results from Ref. ll|. The data are consistent 
with no continuum contribution. 

We report first evidence for T(6S') —>■ /ib(lP)7r+7r“ 
and first observation of the T(6S') —>■ /ib(2P)7r+7r“ tran¬ 
sitions. We study their resonant substructures and find 
evidence that they proceed entirely via the intermediate 
isovector states Zb(10610) and Zf,(10650). Their relative 
fraction is loosely constrained by the current data: the 
hypothesis that only Zb(10610) is produced is excluded at 
the 3.3(7 level, while the hypothesis that only Zb(10650) 
is produced is not excluded at a significant level. 

The shapes of the e“''e“ /ib(nP)7r+7r“ and e+e” ^ 

T(nS')7r'*'7r“ cross sections look similar. The only signif¬ 
icant difference is a smaller relative yield of T(nS')7r+7r“ 
at the T(6S'). Since the /if,(nP)7r“*'7r“ final states are pro¬ 
duced only via intermediate Zb while T(n5')7r“'"7r“ at the 
T(5iS') are produced both via Zb and non-resonantly, this 
difference indicates that the non-resonant contributions 
in T(n5')7r+7r“ are suppressed at the T(6S'). 

We thank the KEKB group for excellent operation 
of the accelerator; the KEK cryogenics group for effi¬ 
cient solenoid operations; and the KEK computer group, 
the Nil, and PNNL/EMSL for valuable computing and 
SINET4 network support. We acknowledge support from 
MEXT, JSPS and Nagoya’s TLPRC (Japan); ARC (Aus¬ 
tralia); FWF (Austria); NSFC and CCEPP (China); 
MSMT (Czechia); CZF, DFG, EXC153, and VS (Ger¬ 
many); DST (India); INFN (Italy); MOE, MSIP, NRF, 
BK21P1US, WCU and RSRI (Korea); MNiSW and NCN 
(Poland); MES, RFAAE and RSF under Grant No. 15- 
12-30014 (Russia); ARRS (Slovenia); IKERBASQUE 
and UPV/EHU (Spain); SNSF (Switzerland); MOE and 
MOST (Taiwan); and DOE and NSF (USA). 


























7 


[1] N. Brambilla et al., Eur. Phys. J. C 71, 1534 (2011). 

[2] S.K. Choi et al. [Belle Collaboration], Phys. Rev. Lett. 
91, 262001 (2003). 

[3] K.A. Olive et al. [Particle Data Group], Chin. Phys. C 
38, 090001 (2014). 

[4] J.P. Lees et al. [BaBar Collaboration], Phys. Rev. D 86, 
051102 (2012); Z.Q. Liu et al. [Belle Collaboration], Phys. 
Rev. Lett. 110, 252002 (2013). 

[5] J.P. Lees et al. [BaBar Collaboration], Phys. Rev. D 89, 
111103 (2014); X. L. Wang et al. [Belle Collaboration], 
Phys. Rev. D 91, 112007 (2015). 

[6] X. L. Wang et al. [Belle Collaboration], Phys. Rev. D 87, 
051101 (2013); M. Ablikim et al. ]BESIII Collaboration], 
Phys. Rev. D 91, 112005 (2015). 

[7] M. Ablikim et al. [BESIII Collaboration], Phys. Rev. 
Lett. Ill, 242001 (2013); C.Z. Yuan, Chin. Phys. C 38, 
043001 (2014). 

[8] N. Brambilla et al, Eur. Phys. J. C 74, 2981 (2014). 

[9] S. Dubynskiy and M. B. Voloshin, Phys. Lett. B 666, 
344 (2008). 

]10] K.-F. Chen et al. ]Belle Collaboration], Phys. Rev. D 82, 
091106 (2010). 

]lll D. Santel et al. ]Belle Collaboration], Phys. Rev. D 93, 
011101 (2016). 

[12] M. B. Voloshin, Phys. Rev. D 85, 034024 (2012). 

[13] 1. Adachi et al. [Belle Collaboration], Phys. Rev. Lett. 


108, 032001 (2012). 

[14] A. Bondar et al. [Belle Collaboration], Phys. Rev. Lett. 
108, 122001 (2012). 

[15] A. Bondar et al., Phys. Rev. D 84, 054010 (2011). 

[16] A. Abashian et al. [Belle Collaboration], Nucl. Instrum. 
Meth. A 479, 117 (2002); also see detector section 
in J. Brodzicka et al., Prog. Theor. Exp. Phys. 2012, 
04D001 (2012). 

[17] S. Kurokawa and E. Kikutani, Nucl. Instrum. Meth. A 
499, 1 (2003), and other papers included in this Volume; 
T. Abe et al.. Prog. Theor. Exp. Phys. 2013, 03A001 
(2013) and references therein. 

[18] R. Mizuk et al. [Belle Collaboration], Phys. Rev. Lett. 

109, 232002 (2012). 

[19] K. Abe et al. [Belle Collaboration], Phys. Rev. D 64, 
072001 (2001). 

[20] G.C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 
(1978). 

[21] E.A. Kuraev and V.S. Fadin, Sov. J. Nucl. Phys. 41, 466 
(1985); M. Benayoun, S.I. Eidelman, V.N. Ivanchenko 
and Z.K. Silagadze, Mod. Phys. Lett. A 14, 2605 (1999). 

[22] A. Garmash et al. [Belle Collaboration], Phys. Rev. D 
91, 072003 (2015). 

[23] S. Actis et al, Eur. Phys. J. C 66, 585 (2010). 

[24] See Supplemental Material at [URL will be inserted by 
publisher] for tables of a^[e'^e~ —> hi,{nP)TT'^'K~] {n = 
1 ,2) measurements. 

[25] S.S. Wilks, Ann. Math. Statist. 9, 60 (1938). 

DOI:10.1214/aoms/1177732360. 



arXiv:1508.06562v2 [hep-ex] 9 Jun 2016 


1 


Supplemental Material 


1 


The e“'"e —>■ hi,{nP)'K~^TT (n 


1, 2) Born cross sections for all energy points are presented in Table HI 


2 TABLE I. Center-of-mass energies, integrated luminosities and Born cross sections for all energy points. The first uncertainty 

3 in the energy is uncorrelated, the second is correlated. The three uncertainties in the cross sections are statistical, uncorrelated 

4 systematic and correlated systematic. 


# 

MeV 

Luminosity, fb ^ 

cr(e+e 


h6(lP)7r+7r 


), Pb 

cr(e+e —>■ /ii,( 2 P) 7 r+ 7 i 


), pb 

1 

11022 . 0 + 5:3 ± 1.0 

0.98 

-0.39 

± 

0.85 

± 

0.45 

± 

0.02 

2.21 ±1.05 ±0.51 

± 

0.15 

2 

11017.5 ±4.0 ±1.0 

0.86 

1.87 

± 

0.90 

± 

0.37 

± 

0.11 

2.05 ± 1.091^:2? 

± 

0.14 

3 

11016.41^:6 ± 1.0 

0.77 

1.25 

± 

0.96 

± 

0.19 

± 

0.08 

3.42± 1.15l°:l3 

± 

0.23 

4 

11006 . 813:9 ± 1.0 

0.98 

2.57 

± 

0.88 

± 

0.13 

± 

0.15 

3.45 ± 1.03 ±0.17 

± 

0.23 

5 

10991.9 ±0.4 ±1.0 

0.99 

2.25 

± 

0.91 

± 

0.16 

± 

0.14 

2.45 ± 1.11 ±0.13 

± 

0.17 

6 

10977.5 ±0.4 ±1.0 

1.00 

0.38 

± 

0.90 

± 

0.33 

± 

0.02 

0.83 ± 1.23 ±0.29 

± 

0.06 

7 

10957.5 ±4.0 ±1.0 

0.97 

1.05 

± 

0.89 

± 

0.14 

± 

0.06 

0.60 ± 1.39lg:l2 

± 

0.04 

8 

10927.5 ±4.0 ±1.0 

1.15 

-0.16 

± 

0.89 

± 

0.13 

± 

0.01 

2.92 ± 1.59 ±0.33 

± 

0.20 

9 

10907.71^:9 ± 1-0 

0.98 

0.97 

± 

0.96 

± 

0.27 

± 

0.06 

-0.04 ± 1.72 ±0.33 

± 

0.00 

10 

10901. 1 II 9 ± 1-0 

1.42 

1.72 

± 

0.79 

± 

0.16 

± 

0.10 

2.76 ± 1.44+0,2° 

± 

0.19 

11 

10898.5l^:o ± 1-0 

0.98 

2.63 

± 

0.96 

± 

0.27 

± 

0.16 

3.89 ± 1.70t^:?^ 

± 

0.26 

12 

10888.9l^:o ± 1-0 

0.99 

1.94 

± 

1.00 

± 

0.48 

± 

0.12 

5.41 ±1.83 ±0.40 

± 

0.37 

13 

10883.61°:? ± 1-0 

1.85 

3.89 

± 

0.76 

± 

0.19 

± 

0.23 

5.55 ± 1.37 ±0.45 

± 

0.38 

14 

10878.5l?:l ± 1-0 

0.98 

2.89 

± 

1.04 

± 

0.23 

± 

0.17 

5.34 ± 1.89 ±0.91 

± 

0.36 

15 

10869.5l?:o ± 1-0 

0.98 

2.22 

± 

1.04 

± 

0.18 

± 

0.13 

4.26 ± 1.90 ±0.68 

± 

0.29 

16 

10868.6 ±0.2 ±0.5 

22.94 

1.36 

± 

0.21 

± 

0.04 

± 

0.08 

3.01 ±0.39 ±0.08 

± 

0.20 

17 

10866.7 ±0.2 ±0.5 

50.47 

1.81 

± 

0.15 

± 

0.04 

± 

0.11 

2.91 ±0.26 ±0.08 

± 

0.20 

18 

10863.3 ±0.2 ±0.5 

47.65 

1.66 

± 

0.15 

± 

0.03 

± 

0.10 

2.33 ±0.27 ±0.11 

± 

0.16 

19 

10858.9l?:o ± 1-0 

0.99 

0.44 

± 

1.03 

± 

0.11 

± 

0.03 

3.63 ± 2 . 011 ^: 2 ? 

± 

0.25 

20 

10849.7l?:2 ± 1.0 

0.99 

1.54 

± 

1.02 

± 

0.13 

± 

0.09 

2.57±2.07l°:S 

± 

0.17 

21 

10820.51^:? ± 1.0 

1.70 

0.50 

± 

0.77 

± 

0.10 

± 

0.03 

1.83 ± 1.40 ±0.70 

± 

0.12 

22 

10771.1 ± 1.8 ± 1.0 

0.95 

0.28 

± 

0.86 

± 

0.19 

± 

0.02 

0.47 ± 1.52 ±0.12 

± 

0.03 





